home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / dlamch.z / dlamch
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. DDDDLLLLAAAAMMMMCCCCHHHH((((3333FFFF))))                                                          DDDDLLLLAAAAMMMMCCCCHHHH((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DLAMCH - determine double precision machine parameters
  10.  
  11. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  12.      DOUBLE PRECISION FUNCTION DLAMCH( CMACH )
  13.  
  14.          CHARACTER    CMACH
  15.  
  16. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  17.      DLAMCH determines double precision machine parameters.
  18.  
  19.  
  20. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  21.      CMACH   (input) CHARACTER*1
  22.              Specifies the value to be returned by DLAMCH:
  23.              = 'E' or 'e',   DLAMCH := eps
  24.              = 'S' or 's ,   DLAMCH := sfmin
  25.              = 'B' or 'b',   DLAMCH := base
  26.              = 'P' or 'p',   DLAMCH := eps*base
  27.              = 'N' or 'n',   DLAMCH := t
  28.              = 'R' or 'r',   DLAMCH := rnd
  29.              = 'M' or 'm',   DLAMCH := emin
  30.              = 'U' or 'u',   DLAMCH := rmin
  31.              = 'L' or 'l',   DLAMCH := emax
  32.              = 'O' or 'o',   DLAMCH := rmax
  33.  
  34.              where
  35.  
  36.      eps   = relative machine precision
  37.            sfmin = safe minimum, such that 1/sfmin does not overflow base  =
  38.            base of the machine prec  = eps*base t     = number of (base)
  39.            digits in the mantissa rnd   = 1.0 when rounding occurs in
  40.            addition, 0.0 otherwise emin  = minimum exponent before (gradual)
  41.            underflow rmin  = underflow threshold - base**(emin-1) emax  =
  42.            largest exponent before overflow rmax  = overflow threshold  -
  43.            (base**emax)*(1-eps)
  44.  
  45.  
  46.  
  47.  
  48.  
  49.  
  50.  
  51.  
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.